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This  report  documents  the  last  of  four  tidal  perturbing  terms  appearing  in  the 
force  equations  for  satellite  motion  in  the  TERRA  system  of  satellite  geodesy 
computer  programs.  An  algorithm  is  presented  for  the  gravitational  action  by  which 
the  tidal  bulge  of  the  solid  earth  perturbs  satellite  motion.  Three  earlier  reports 
contain  similar  algorithms  for  the  air  tides  caused  by  the  sun  and  the  moon,  as  well 
as  for  the  ocean  tides. 

The  problem  analysis  was  rather  tedious  because  it  involved  certain  expressions 
which  were  unusually  complicated  algebraically.  The  authors  are  greatly  obliged  to 
Mr.  Alfred  Morris  of  the  Computer  Program  Division  (DK-70),  who  verified  the 
results  with  the  aid  of  his  symbolic  computer  language  FLAP. 

The  work  was  done  in  the  Astronautics  and  Geodesy  Division  (DK-10)  in 
support  of  astronautics  computer  program  development.  This  report  has  been 
reviewed  and  approved  by  R.  J.  Anderle,  Head,  Astronautics  and  Geodesy  Division. 
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INTRODUCTION 


All  tides  except  one’  are  caused  by  the  fields  of  lunar  and  solar  mass 

attraction  within  which  the  earth  rotates.  These  fields  impart  to  each  terrestrial 

particle  an  acceleration  which  depends  on  the  specific  position  of  the  particle  within 
the  field.  Generally,  this  acceleration  is  different  from  that  experienced  by  the  earth 
as  a whole.  The  relative  acceleration  of  the  particle  with  respect  to  the  earth  gives 
rise  to  certain  geophysical  forces  which  balance  it,  thus  enabling  the  particle  to 

remain  attached  to  the  earth.  These  forces  constitute  a very  complicated  stress  field. 
It  is  customary  to  classify  the  various  tides  by  decomposing  the  tidal  stress  field 

into  a number  of  individual  components  according  to  the  nature  of  the  terrestrial 
constituents  (solid  earth,  ocean,  or  atmosphere)  and  to  the  body  (sun  or  moon) 
which  generates  the  particular  stress  field  component. 

Each  tide  is  essentially  a redistribution  of  the  earth’s  masses.  As  such,  it  will 
modify  the  terrestrial  gravity  field  which  in  turn  may  be  expected  to  perturb  the 

satellite  motion.  For  most  purposes  of  satellite  geodesy,  tidal  perturbation  of  the 
orbits  was  regarded  inconsequential  until  quite  recently.  However,  data  acquisition 
has  now  become  so  accurate  and  the  art  of  orbit  computation  has  been  refined  to 
the  point  where  it  is  generally  agreed  that  to  avoid  tide  related  biases,  the  tides  are 
required  to  appear  amongst  the  perturbing  terms  in  the  force  equations  for  satellite 
motion. 

The  necessary  improvement  in  the  force  model  for  geodetic  satellites  was 
recently  implemented  at  the  Naval  Surface  Weapons  Center.  The  additional  force 
terms  which  have  been  developed  for  the  new  TERRA  system  of  computer  programs 
for  satellite  geodesy  reflect  the  state  of  the  art  and  account  for  tides  due  to  the 
combined  effects  of  solar  gravitation  and  solar  irradation,1  tides  caused  by  the 

moon’s  gravity,  ocean  tides,3  and  solar  as  well  as  lunar  tides  of  the  solid  earth. 
For  the  two  atmospheric  tides,  surface  pressure  functions  corresponding  to  the 

atmospheric  tide  bulges  were  readily  available  in  the  geophysical  literature.  The 
disturbing  acceleration  was  derived  for  each  tide  from  surface  pressure  via  Poisson 

integration  followed  by  spatial  differentiation.  The  algorithm  for  the  ocean  tides 
starts  from  a table  for  the  global  tide  amplitudes  and  phase  angles.  It  approximates 
the  gravitational  tide  potential  exterior  to  the  earth  by  that  of  a system  of  point 
masses  which  are  located  on  the  ocean  surface  and  which  are  chosen  to  resemble 
the  masses  of  the  tidal  bulge. 


* fhc  solar  air  tide  i>  partly  a thermal  effect  related  to  the  absorption  ol  sunlight. 
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A suitable  potential  function  for  the  solid  earth  tide  was  compiled  from  current 
literature.4  This  was,  however,  very  complicated  algebraically.  The  partial  derivatives 
necessary  for  the  perturbing  acceleration  were  derived  manually;  and  to  assure 
correct  results,  this  effort  was  duplicated  on  the  computer  with  the  help  of  the 
symbolic  language  FLAP.5,6  In  addition,  it  was  difficult  to  verify  by  hand,  in  a 
reasonable  time  span,  that  the  formula  collection  for  the  potential  actually  satisfied 
the  Laplace  differential  equation.  Consequently,  verification  was  performed  entirely 
by  computer  which  speedily  resolved  the  problem. 

The  potential  function  just  mentioned  is  probably  the  most  flexible  expression 
now  available  for  the  Newtonian  gravitation  associated  with  the  solid  earth  tide. 
However,  it  was  not  immediately  applicable  to  our  purpose  because,  as  transcribed 
from  Reference  4,  it  represents  the  mass  attraction  of  the  equilibrium  tide.  Our 
attention  was  directed  to  this  defect  by  Reference  7 (page  2)  which  also  indicates 
that  the  necessary  correction  can  be  accomplished  by  interpreting  the  just  mentioned 
equilibrium  tide  potential  as  being  produced  by  a fictitious  moon  or  sun.  A rather 
simple  transformation  from  the  lunar  and  solar  ephemeris  coordinates  to  those  of 
the  fictitious  moon  or  sun  will  change  the  equilibrium  potential  into  a potential 
function  which  contains  the  desired  effects  of  tidal  lag. 

THE  TIDE  POTENTIAL 


The  following  is  a summary  of  the  tide  potential  associated  with  the  tidal  mass 
redistribution  occurring  within  the  solid  earth.  This  was  abstracted  from  Reference  4 
(see  especially  pages  258  through  266).  The  formulas  chosen  were,  however,  not 
simply  copied.  When  comparing  the  equations  which  appear  on  the  next  pages  with 
their  counterparts  in  the  reference,  one  will  notice  that  the  unit  vector  components 
and,  generally,  the  coordinates  which  indicate  respectively  the  positions  of  sun  and 
moon  are  distinguished  below  from  those  of  the  satellite  by  a star-shaped 
superscript,  while  in  the  reference,  they  are  primed.  In  the  reference  as  well  as  here, 
the  prime  denotes  the  ephemeris  coordinates,  while  the  star  is  attached  to  the 
positions  of  the  fictitious  sun  and  moon  mentioned  above  in  the  Introduction  and 
further  elaborated  in  the  next  chapter. 


*This  it  j stationary  solution  of  the  tidal  problem.  If  taken  for  the  tide  proper,  the  unrealistic  assumption  would  be 
implied  of  a tidal  mass  redistribution  capable  of  instantaneously  responding  to  the  motion  of  the  celestial  body  which 
causes  the  tide. 
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The  potential  is.  of  course,  valid  exterior  to  the  earth.  Observe  that  the 
properties  of  the  tide  are  assumed  to  vary  latitudinally;  that  is,  they  are  expressed 
by  including  Love  coefficients  which  depend  on  latitude.  To  be  precise,  the  Love 
numbers  are  expansions  in  terms  of  zonal  harmonics,  latitude’  being  the  variable. 
Although  the  expansion  coefficients  for  the  Love  numbers  occur  frequently  in  the 
tide  potential  and  later  in  the  disturbing  acceleration,  the  explicit  expressions  for  the 
Love  numbers  are  not  needed.  Those  interested  in  the  geophysical  background  may 
wish  to  peruse  pages  200  and  261  of  Reference  4.  An  introduction  to  the  basic 
Love  number  concept  may  be  found  in  Reference  8. 

In  detail,  the  tide  potential  is: 
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*We  intentionally  avoid  distinguishing  between  the  various  types  of  latitude  because  we  expect  that  any  inaccuracy  arising 
from  this  simplification  will  be  masked  by  the  much  larger  uncertainties  which  otherwise  prevail  in  current  tide  work. 
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The  k(j  terms  occurring  in  the  formulas  above  are  the  expansion  coefficients  for 
the  Love  numbers.  These  are  geophysical  constants  of  the  TERRA  system.  It  is 
expected  that  the  following  values  approximate  the  k(j  sufficiently;  otherwise,  they 
may  be  regarded  to  be  initial  values  for  a parameter  improvement  procedure. 


k20 

= 0.3 

k2  1 

= 0.01 

k30 

= 0.1 

k3  1 

= 0.01 

\ 

. 

(10) 

k2  2 

= 0.1 

i 

The  necessary  notations  are  compiled  in  Table  I. 
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Table  1.  Glossary  for  the  Tide  Potential  Formulas 


Mass  of  the  earth 
Newton’s  constant 

Equatorial  radius  of  the  earth  (semimajor  axis  of  the  reference 
ellipsoid  in  use) 

Eccentricity  of  earth’s  meridian  (eccentricity  of  the  reference  ellipsoid 
in  use) 

Expansion  coefficients  for  the  Love  numbers 
Mass  of  the  sun  or  moon 

Semimajor  axis  of  the  sun  or  moon  orbit  relative  to  the  earth 
(Numerical  values  to  be  compatible  with  the  ephemeris  tapes  used  in 
TERRA  for  the  sun  and  moon.  May  optionally  be  taken  from  a 
current  edition  of  “American  Ephemeris”.) 

Equatorial  components  of  the  unit  vector  ^ pointing  from  the  center 
of  the  earth  to  the  fictitious  sun  or  moon  in  a Cartesian  coordinate 
system  referred  to  mean  equator  and  equinox  of  epoch 

Geocentric  distance  of  fictitious  sun  or  fictitious  moon 

Mean  motion  of  the  satellite 

Semimajor  axis  of  the  satellite  orbit  (This  is  an  osculating  orbit 
element  which  may  be  obtained  from  postion  and  velocity  of  the 
satellite,  for  each  time  line  in  the  orbit  integration,  via  the  energy 
integral.) 

Equatorial  components  of  the  unit  vector  ? pointing  from  the  center 
of  the  earth  to  the  satellite  in  a Cartesian  coordinate  system  referred 
to  mean  equator  and  equinox  of  epoch 


Geocentric  distance  of  the  satellite 


TIDAL  LAG 


We  mentioned  earlier  that  the  original  equations  for  the  tide  potential 
(transcribed  from  Reference  4)  are  part  of  a stationary  tide  model,  bssential  to  the 
latter  is  the  assumption  that  the  tidal  mass  redistribution  occurs  instantaneously  or 
(which  is  the  same)  that  the  motions  relative  to  the  earth,  of  sun  and  moon, 
involve  but  vanishing  velocity.  Hither  one  of  these  notions  is  unrealistic  because  they 
omit  tidal  lag.  lor  THRRA,  we  desired  to  include  the  latter  effect,  hoping  to  be 
able  to  measure  how  it  acts  on  the  satellite  orbits. 

The  tidal  lag  feature  was  readily  introduced  by  the  expedient  of  a fictitious 
sun  and  a fictitious  moon.  These  are  the  ephemeris  sun  and  the  ephemeris  moon: 
each,  however,  retarded  in  their  geocentric  motion  with  respect  to  the  ephemeris 
positions  and  advanced  in  the  direction  of  the  earth’s  rotation.  The  retardation  in 
the  geocentric  motion  is  expected  to  reflect  how,  at  any  given  location  on  the 
globe,  viscosity  and  inertia  combine  to  delay  the  forming  of  the  tidal  bulge  while 

the  sun  and  moon  traverse  the  sky.  Simultaneously,  the  advance  of  either  of  the 

two  fictitious  bodies  in  the  direction  of  the  earth’s  rotation  represents  the  tidal 
bulge  being  carried  along  by  that  rotation. 

As  illustrated  in  Figure  1,  the  motion  geometry  just  specified  was  realized  by  a 
transformation  from  the  ephemeris  coordinates  of  the  sun  and  moon  to  the 
corresponding  fictitious  positions.  Separately,  for  each  of  the  two  bodies,  the  orbital 
retardation  was  accomplished  by  including  a time  delay  (equal  to  the  tidal  lag 

parameter  appearing  in  the  formulas  below)  into  the  interpolations  which  produce 
the  instantaneous  sun  and  moon  positions  from  the  ephemeris  data.  By  contrast,  the 
rotational  advance  was  introduced  by  a suitable  rotation  of  the  orbital  planes 
associated  with  the  motions,  relative  to  the  earth,  of  the  sun  and  moon.  The 

following  algorithm  resulted  for  the  combined  effect  of  the  two  motions:  and,  the 
notation  is  explained  in  Table  2. 
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Table  2.  Glossary  for  the  Tidal  Lag  Equations 


V yr  zj  Table  values  of  the  coordinates  for  the  sun  and  moon  from  the 

ephemeris  tapes 

•AT  Time  increment  from  the  ephemeris  tapes 

At  Tidal  lag  parameter 

x\  y\  z Position  of  the  sun  or  moon  at  time  t as  computed  from  the 

ephemeris  tape  by  interpolation  (“ephemeris  coordinates”) 

x*,  y*,  z*  Coordinates  of  the  fictitious  sun  or  moon  at  time  t 

n#  Earth’s  sidereal  rate  of  rotation 


i 


* UNPI  RUIRHt  0 ORBl!  III  SUN  MIH1S 
B PIBIUHBID  ORBIT  01  SUN  MOON 
P [ PHI  Ml  HIS  POSITION  01  SUN  MOON 
P'  Ri  TAROtO  POSITION  01  SUN  MOON 


Figure  1.  Simulating  Tidal  Lag  by  Orbital  Retardation 
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TIDES  OF  THE  SOLID  EARTH  IN  THE  TERRA  EQUATIONS  OF  MOTION 


The  desired  perturbing  term  can  now  readily  be  formulated.  Note  that  the 
vector  T just  introduced  (Equation  16)  depends  essentially  on  the  coordinates  ol  the 
fictitious  sun  or  moon,  on  the  mass  and  the  semimajor  axis  of  the  sun  or  moon, 
and  on  the  satellite  coordinates: 


T = T(x*,y*,  z*,m\a'  jx^y^z.) 


( 18) 


To  produce  the  perturbing  acceleration  caused  by  the  solar  tide,  substitute  into 
the  latter  equation,  separately  for  each  step  (time  line)  in  the  satellite  orbit 
integration,  the  position  coordinates  of  the  fictitious  sun,  the  satellite  position,  and 
the  remaining  two  solar  parameters  (mj.  and  a^ ): 


Ts  = Ts<xs-ys’zs*’ms'as  ; X,,  y.,  Zj) 


(19) 


For  the  lunar  tide  component,  enter  the  corresponding  (lunar)  data  to  obtain: 


TL  = VXlV^’mL’aL  ; xj.  Vj.  Zj> 


(20) 


The  tide  potential,  as  compiled  by  Reference  4,  is  the  sum  of  the  potentials 
associated  with  the  lunar  and  solar  tide  components.  This  implies  that  the  disturbing 
acceleration  due  to  the  combined  effect  of  the  two  tide  components  is,  simply,  the 
sum  of  the  individual  accelerations: 


’SOI.  10  KART  II 


= T 

t i o k:  s 


+ T, 


(21) 


All  the  equations  above  are  referred  to  an  inertial  reference  frame.  The  latter  is 
identical  with  the  coordinate  system  in  terms  of  which  TERRA  performs  the 
satellite  orbit  integration.  In  particular,  the  perturbing  acceleration  specified  by 
Equation  21  is  compatible  with  all  the  other  perturbing  terms  appearing  in  the 
TERRA  equations  of  motion.  No  further  coordinate  transformation  is  needed.  The 
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force  equations  for  satellite  motion  in  TERRA  are  a rather  complicated  collection  of 
algorithms,  each  of  which  reflects  one  of  the  perturbing  forces  acting  on  earth 
satellites.  The  entire  force  model  has  been  formulated  and  is  available  for 
implementation  on  the  computer;  and,  it  will  be  discussed  in  a future  report. 

It  should  finally  be  noted  that  the  tide  potential,  as  it  appears  in  Reference  4, 
is  a function  of  the  ephemeris  coordinates  of  moon  and  sun.  To  introduce  tidal  lag, 
we  substituted  for  the  ephemeris  positions  the  fictitious  positions  of  moon  and  sun. 
We  are  aware  that  this  is  not  strictly  correct.  To  insert  the  tidal  lag  feature  in  the 
right  place,  it  would  have  been  necessary  to  re-derive  the  tide  potential.  The 
positions  of  the  fictitious  moon  and  sun  could  then  have  been  considered  while 

dealing  with  the  tidal  surface  deformation,  in  the  earth-fixed  reference  frame,  prior 
to  calculating  the  tide  potential  and  prior  to  transforming  the  latter  back  into  the 
inertial  frame.  This  would  have  accurately  accounted  for  the  effects  of  precession 
and  nutation.  But  as  we  expect  these  to  be  negligible  when  compared  with  the  lack 
ol  precision  inherent  in  our  tide  lag  model,  we  decided  on  the  approach  actually 
chosen. 

The  reader  may  now  expect  to  be  presented  with  a quantity  of  numerical 

information  concerning  items  such  as  the  average  magnitude  of  the  perturbing 
acceleration  at  typical  orbital  altitude,  plots  of  the  disturbing  force  along  a typical 
satellite  trajectory,  and  estimates  of  the  tidal  effect  on  the  various  orbital  elements. 
However,  no  serious  attempt  was  made  to  produce  such  data  by  manual  calculation. 

Rather,  the  authors  are  planning  to  take  advantage  of  the  fact  that  the  TKRRA 

force  equations  are  intended  for  use  with  a direct,  numerical,  orbit  integration.  It  is 
expected  that,  upon  completion  of  the  program  coding,  all  numerical  questions  will 
best  be  explored  by  trial  runs  of  the  relevant  program  segments. 
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